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ABSTRACT 

In  this  work,  the  high-order  central  compact  scheme  is  introduced  to 
improve  the  accuracy  of  the  IFLOW  solver.  The  high-order  compact  filter 
is  used  to  reduce  non-physical  oscillations.  To  achieve  the  conservation 
property  of  the  scheme,  the  compact  scheme  is  combined  with  the  ENO 
reconstruction  method.  The  numerical  flux  at  the  cell  interface  is 
approximated  by  the  derivative  of  its  primitive  function  which  can  be 
exactly  calculated.  At  the  boundary,  ghost  points  are  constructed.  The 
high-order  extrapolation  is  used  to  obtain  values  at  the  ghost  points. 
Instead  of  using  the  one-sided  compact  scheme,  the  explicit  central 
difference  scheme  is  applied  at  the  boundary.  The  results  from  spatial 
simulations  of  trailing  vortices  have  shown  obvious  improvement  in 
accuracy. 

1.  INTRODUCTION 

The  compact  schemes  have  been  widely  used  in  numerical 
simulation  of  complex  flows.  It  has  been  shown  [1]  that  the  compact 
scheme  can  reach  high-order  accuracy  and  good  small  scale  resolution 
with  narrow  grid  stencils. 

IFLOW  code  [2]  is  designed  as  a general-purpose  production  code 
for  the  numerical  solutions  of  the  incompressible  Reynolds-averaged 
Navier-Stokes  equations  supplemented  by  appropriate  turbulence  models. 
The  numerical  methods  include  multi-block  grid  structure  for  complex 
geometry,  artificial  compressibility  model,  explicit  one-step  multi-stage 
Runge-Kutta  scheme  for  time  stepping,  mutigrid,  local  time-stepping, 
implicit  residual  smoothing,  optimal  precondition,  and  bulk  viscosity 
damping  for  convergence  acceleration.  The  second  order  central 
difference  and  finite  volume  method  are  used  for  spatial  discretization. 
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The  fourth-order  artifical  dissipation  terms  are  added  to  suppress 
numerical  oscillations.  The  code  has  been  used  to  analyze  many  complex 
flow  fields  including  those  around  submarine,  surface  ship,  and  various 
parts  of  each. 

In  this  work,  the  fourth-order  compact  finite  difference  scheme  is 
incorporated  into  the  code  to  improve  the  accuracy  of  the  solutions. 
Instead  of  using  artificial  dissipation,  the  sixth  order  compact  filter  is  used 
to  reduce  non-physical  oscillations.  To  achieve  the  conservation  property 
of  the  scheme,  the  compact  scheme  is  combined  with  the  ENO 
reconstruction  method[3].  The  numerical  flux  at  the  cell  interface  is 
approximated  by  the  derivative  of  its  primitive  function  which  can  be 
exactly  calculated.  The  physical  flux  at  the  boundary  is  used  as  the 
boundary  condition.  To  calculate  the  physical  flux,  ghost  points  are 
constructed  by  the  fourth  order  extrapolation.  The  next  section  presents 
this  method  in  detail,  followed  by  the  test  cases  and  results  of  calculations. 


2.  GOVERNING  EQUATIONS 


The  three-dimensional  incompressible  Reynolds-averaged  Navier- 
Stokes  equations  based  on  the  artificial  compressibility  approach  can  be 
expressed  using  the  following  conservative  formulation. 


Po'q, 


+ Fx  +Gy  + H z = 0 


(1) 


where  the  preconditioned  matrix  P0  and  the  three  components  of  fluxes 
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Here  p*  is  the  pressure  p divided  by  a constant  density  p . a , ft  2 and 
y are  preconditioning  parameters.  In  this  presented  work,  a = y = 0 and 
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/3‘2  = max(|  w|\  £ ) £ = 0.7.  The  Reynolds  stresses  calculated  by 
turbulence  model  are  given  by  t0  , i,  j = x,y,z- 


3.  NUMERICAL  METHOD 


In  the  old  version  of  IFLOW,  a finite  volume  method  is  used.  The 
mean  flow  is  descretized  by  a second-order  accurate  central  difference 
method  with  the  fourth-order  dissipation  terms.  The  turbulent  flow  is 
discretized  by  one  of  the  several  upwind  schemes.  The  time  step  is  based 
on  an  explicit  one-step  multi-stage  Runge-Kutta  method  to  reach  a steady- 
state  solution.  Several  convergence  acceleration  techniques  including 
multigrid,  local  time  step,  implicit  residual  smoothing,  preconditioning, 
and  bulk  viscosity  damping  have  been  implemented. 

In  this  work,  the  fourth-order  compact  finite  difference  scheme  [1] 
and  the  ENO  reconstruction  method  [3]  are  used  to  approximate  the 
derivatives  of  fluxes.  To  illustrate  how  to  use  this  method,  we  choose  Fx , 
the  derivative  of  flux  in  the  * direction,  as  the  example.  When  a 
conservative  approximation  to  the  spatial  derivative  is  applied,  Fx  can  be 


expressed  as 


_1_ 

Ax 


( 


\ 

) 


(2) 


where  F { and  F y are  numerical  flux  functions  at  the  cell  interfaces.  In 

;+2  J~2 

order  to  achieve  the  high-order  accuracy,  numerical  fluxes  should  be 
defined  in  such  a way  that  the  difference  of  numerical  fluxes  is  a high- 
order  approximation  of  the  derivative  Fx . According  to  the  ENO 
reconstruction  procedure  [3],  it  has  been  proved  that  the  primitive  function 
of  F at  the  cell  interfaces  can  be  exactly  calculated  by  the  value  of  F . at 


given  points.  If  P is  the  primitive  function  of  F , then: 
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Then  the  numerical  flux  F at  the  cell  interface  can  be  obtained  by 
taking  derivative  of  its  primitive  function  P, 

f , = p; , (4) 
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The  approximations  of  derivatives  of  P at  cell  interfaces  are  calculated  by 
the  fourth-order  compact  finite  difference  scheme.  Gy  and  H z are 

discretized  in  the  same  manner.  The  sixth  order  compact  filter  is  used  to 
reduce  numerical  oscillations. 

4.NUMERICAL  EXAMPLES 

To  verify  the  efficiency  of  this  high-order  approach,  the 
development  of  a trailing  vortex  is  simulated  using  the  k-ca  turbulence 
model.  In  the  first  testing  problem,  a laminar  q-vortex  solution  of 
Batchelor[2]  is  used  as  the  inlet  flow  condition,  i.e. 

1 rT  ( r2YI  ( r2} 

+ — )-L  1-exp  -a—  U-U_=U„  -a-j 

20  rY  'i  JJ  d J 

where  r2  = y2  + z2 , a = 1 .25643 , Vai  =0.286,  UD  =0.165,  and 

r,  = 0.036 . The  size  of  the  computational  domain  is  40  x 1 x 1 . The  inlet  is 
located  at  jc=5.  The  grid  number  is  96x32x32.  Fig.l  and  2 show  the 
distributions  of  the  streamwise  and  tangential  velocities  along  a vertical 
line  through  the  vortex  core  at  different  locations  in  the  x direction.  The 
results  from  the  second-order  and  the  fourth-order  are  displayed  in  the 
same  figure.  Fig.3  shows  the  contours  of  vorticity  at  different  x locations. 
The  first  row  is  corresponding  to  the  second-order  results,  and  the  second 
row  is  of  the  fourth-order.  It  can  be  seen  from  these  results  that  the  high- 
order  compact  scheme  increases  the  accuracy  of  the  solutions  and  has 
demonstrated  a better  resolution  than  second-order  scheme. 


Fig.l  Streamwise  velocity  Fig.2  Tangential  velocity 
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Fig.3  Contours  of  vorticity  at  different  x locations 


In  the  second  testing  problem,  the  experimental  data  from  [2]  are  used  to 
prescribe  the  inlet  flow.  The  geometry  and  flow  condition  match  those  in 
the  experiment.  The  computational  domain  extends  from  x/c=5  to  x/c=45 
in  the  streamwise  direction.  The  grid  is  clustered  near  the  center  of  the 
vortex  on  the  cross  section.  In  the  streamwise  direction  the  grid  is 
uniform.  The  grid  number  is  56  x 56  x 56 . Experimental  data  ( 17  x 17 ) at 
x/c=5  are  interpolated  to  the  computational  grid  nodes.  The  numerical 
results  are  compared  with  experimental  data  at  x/c=10.  Fig.4  and  5 show 
the  distributions  of  the  streamwise  and  tangential  velocities  along  a 
vertical  line  through  the  vortex  core  at  different  locations  in  the  x 
direction.  The  results  from  the  second-order  and  the  fourth-order  are 
compared  with  the  experimental  results.  It  is  obvious  that  the  accuracy  is 
improved  by  using  the  high-order  compact  scheme.  The  large  discrepancy 
between  numerical  results  and  experiment  data  is  caused  by  errors  of  the 
interpolated  inlet  boundary  condition,  as  we  find  this  discrepancy  also 
exists  at  the  sections  very  close  to  the  inlet  boundary. 

5.CONCLUSIONS 

The  fourth-order  compact  scheme  in  the  conservative  form  has  been 
incorporated  in  to  the  production  code  IFLOW.  The  results  from  the 
testing  cases  have  shown  the  improvement  in  accuracy.  The  further 
applications  of  this  method  to  more  complex  flows  are  undertaken. 
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Fig.4  Streamwise  velocity  Fig.5  Tangential  velocity 
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